Quantum Oscillations in a Topological Insulator Bi 2 Te 2 Se 
with Large Bulk Resistivity (6 ficm) 



O 

(N 

G 



i 

C 

o 
o 



> 
in 

m 



x 



Jun Xiong 1 , A. C. Petersen 1,2 , Dongxia Qu 1 , R. J. Cava 2 and N. P. Ong 1 

1 Department of Physics and 2 Department of Chemistry, Princeton University, New Jersey 08544, U.S.A. 

(Dated: January 10, 2011) 

We report the observation of prominent Shubnikov-de Haas oscillations in a Topological Insulator, 
Bi2Te2Se, with large bulk resistivity (6 f2cm at 4 K). By fitting the SdH oscillations, we infer a large 
metallicity parameter kpi = 41, with a surface mobility (p 3 ~ 2,800 cm 2 /Vs) much larger than the 
bulk mobility fjiti, ~ 50 cm 2 /Vs). The plot of the index fields B v vs. filling factor v shows a |-shift, 
consistent with massless, Dirac states. Evidence for fractional-filling states is seen in an 11-T field. 



Topological Insulators are predicted to bear current- 
carrying, massless, Dirac surface states that traverse 
the bulk energy gap jl'j-4]. These unusual surface 
states have been observed by angle-resolved photoemis- 
sion spectroscopy (ARPES) 0tM| and scanning tunnel- 
ing microscopy experiments (STM) [9]. Quantization 
of the Dirac states into Landau Levels has been demon- 
strated in STM experiments [13, [HI ■ Observation of the 
surface currents by transport has been more challeng- 
ing [12) Hit. Recently, however, encouraging progress has 
been achieved. The surface SdH oscillations and surface 
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mobility was measured in Bi2Te3 crystals 
istence of surface states at fractional filling in a pulsed 
magnetic field were reported in (Bi,Sb)Se3 [15] . Here 
we report the observation of prominent SdH oscillations 
in crystals of Bi2Te2Se with very large bulk resisitivity. 
Evidence for fractional-filling states become apparent at 
relatively low magnetic fields. 

Recently, ARPES experiments have shown that 
Bi2Te2Se displays a single topological surface state [l6j . 
Independent of our experiment, SdH oscillations in 
this compound have also been reported by Y. Ando's 
group 171. We were motivated to grow crystals of 
Bi 2 Te2Se by reasoning that Bi 2 Se3 crystals grow with 
a stable density of Se vacancies. Electrons donated by 
the vacancies pin the Fermi energy Ep to the conduction 
band, resulting in a negative thermopower S at low tem- 
perature T. By contrast, as-grown Bi2Te3 typically has 
a positive thermopower. By growing a series of the hy- 
brid semiconductor Bi2Sei+ 2; Te2- :I :, we have found that 
the low-temperature thermopower varies systematically, 
reflecting changes in Ep. Crystals were grown by a modi- 
fied Bridgeman method from high purity elemental start- 
ing materials. After heating for one day at 850 C in a 
clean evacuated quartz tube, the melt was cooled in a 
temperature gradient to 500 C where it was left to an- 
neal for 2 days before cooling rapidly to room tempera- 
ture. For the composition Bi2Te2Se, S rises to very large 
values. With a razor blade, we cleaved thin crystals and 
attached contacts using silver paint. For the sample here, 
the crystal thickness d = 110 /im, while the distance be- 
tween voltage leads equals 0.5 mm. The resistivity profile 
measured at B = is shown in Fig. [1^. Below 40 K, the 
value of p attains values in the range 5-6 51cm, or ~1000 
times higher than in non- metallic Bi2Te3. Expressed as 
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FIG. 1: (Color online) The resistivity of a cleaved crystal of 
the Topological Insulator Bi2Te2Se. Panel (a) shows p vs. T 
measured in B = 0. Below 10 K, p attains the value 5.5 ficm, 
or an areal resistance Rq = 400 Q. Despite the non-metallic 
value of Ra, sizeable quantum oscillations are observed be- 
low 38 K. Panel (b) displays the prominent SdH oscillations 
observed in the derivative dp xy /dB vs. B at 4.4 K. 



an areal resistance iZp = p/d, the low-T resistance cor- 
responds to i?n = 400 f2. Despite the high resistance, p 
is only weakly T-dependent, displaying a logT increase 
as T — > 0. The observed Hall coefficient Rjj below 10 K 
(n-type) implies a very small bulk carrier concentration 
rib ~ 2.6 x 10 16 cm -3 . Combining this with the observed 
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p, we infer a low bulk mobility pb ~ 50 cm 2 /Vs. 

Such a low Hb should not produce SdH oscillations for 
B <14 T. Surprisin gly, however, the Hall resistivity p yx 
displays prominent SdH oscillations that may be resolved 
up to 40 K. To date, we have detected SdH oscillations 
in 4 crystals of Bi2Te2Se with p-T profiles similar to that 
in Fig. QJi. Figure QJ> shows the trace of the derivative 
dp yx /dT vs. B at T = 4.4 K. Independently, SdH os- 
cillations were also observed in Bi2Te2Se by Y. Ando et 
al. [13 ■ 

The surface conductance and bulk conductance act as 
parallel channels for charge transport. The observed con- 
ductivity dij is then the sum 
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where is the bulk conductivity and G?|- the conduc- 
tance matrix of the surface states. To exploit the addi- 
tivity, we have converted the measured resistivity matrix 
Pij to the conductivity matrix er^. To amplify the sur- 
face contribution to a xv , we define Aa xy = a xy — {(J xy ) 
where (a xy ) is a smoothed background. 

Figure [2] displays the subtracted quantity Aa xy ver- 
sus 1/B over the temperature interval 0.3 < T < 38 K. 
With increasing T, the amplitude of the SdH oscillations 
decreases rapidly as T increases above 10 K. 

To analyze the SdH oscillations, we have fitted the 
curves using the standard expression [l8[ 
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with A = 27r 2 fcsT/fio; c and \d = 2% IzbTu /huj c , where 
w c is the cyclotron frequency and the Dingle tempera- 
ture is given by To = h/(2nkBT), with r the lifetime. 
Compared with the SdH expression for the conductiv- 
ity <Jxx, the phase </> in the Hall conductivity is shifted 
by 7r/2 (— 7r/4 — > 7r/4). For 2D systems, we may write 
the SdH frequency as 2irEpB / (?uu c ), which simplifies 
to 4ir 2 hn s /e, with the 2D carrier density n s — k 2 F /An 
(per spin) . As shown in Ref. , Eq. [5] may be em- 
ployed in a Dirac system if we write the cyclotron mass 
as m c = E/vp. 



We have found that the oscillations cannot be fitted 
using one SdH frequency. For T < 6 K, the sharp de- 
crease of the oscillation amplitude for B^ 1 > 0.12 T _1 
suggests beating between 2 terms of nearly equal frequen- 
cies. Indeed, a good fit is obtained if we add a second 
term identical to the one in Eq. [2] except for a slight 
difference in n s . The measured curve of Aa xy was fitted 
with the 2 terms (the absolute value of the surface Hall 
conductance G xy is not known). There are altogether 5 
adjustable parameters (n s i and amplitude Ai, with i 
1,2) and To (assumed same for both). The best fit (Fig. 
|3^,) is obtained with A\ — and and densities differing 
by only 5% [(n sl ,n s2 ) = (1.79, 1.71) xlO 12 cur 2 ], cor- 
responding to an average Fermi wavevector kp = 0.047 
A -1 . The fit yields Tp> = 8. 5 ±1.5 K, which corresponds 



FIG. 2: (Color online) The Hall SdH oscillations versus 1/B at 
selected T. In Panel (a), we have plotted the subtracted Hall 
conductivity Aa xy to highlight the decrease of the amplitude 
over 2 decades in T (0.3 to 38 K). Panel (b) displays traces of 
the Hall resistivity dp yx /dB at 5 temperatures. The minima 
in \dp yx /dB\ are used to fix the index field B v (dashed lines 
with v indicated). For LL with v <6, sharp structures are 
observed between integer filling. 



to a mean-free-path £ = 70-100 nm and a surface mobil- 
ity p s — el/hkp — 2,800±250 cm 2 /Vs. Interestingly, p s 
is strongly enhanced over pi, (by ~60). 

By fitting to the decrease in amplitude with T at fixed 
B (12 .4 T), we obtain an effective mass m* = 0.089 m e . 
Together with kp, we obtain a Fermi velocity vf ~ 6 
xlO 5 m/s, higher than that in Bi2Te3. 

The high mobility provides strong evidence that the 
SdH oscillations arise from surface states. Suppose for 
the sake of argument that the oscillations come from bulk 
states. The SdH period must then be identified with a 3D 
Fermi sphere of radius kp — 0.047 A -1 , or a 3D density 
= of 3.3xl0 18 cm~ 3 . The inferred p then implies a 3D 
resistivity pb ~ 0.7 mSlcm at 4.4 K. The large discrepancy 
(factor of 9,000) from the observed value argues firmly 
against a bulk origin. 

Hence we conclude that the SdH oscillations come 
from high mobility surface carriers. The two periods 
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FIG. 4: (Color online) Construction used to fix the filling 
factor for a Dirac spectrum. Panel (a) sketches schematically 
the step-like increase of a xy — (e 2 /h)[v + |] versus energy 
E, where ^ arises from the v = LL at the Dirac point. 
The half-ovals are broadened Landau Levels centered at E v — 
hig Vfv2v. B v are the fields at which £ falls between LLs 
(arrows) . The inset shows the 2D Dirac energy surface in zero 
B. In Panel (b), the measured values of \jB v are plotted 
versus filling factor v. The 5 values of 1/B„ fall on a straight 
line that intercepts the v axis at -0.55, consistent with a Dirac 
spectrum. 



FIG. 3: (Color online) Fits of Aa xy to extract mobility. Panel 
(a) shows the fit to Aa xy at T = 4.4 K to Eq. The 
rapid decrease of the oscillation amplitude for 1/B > 0.09 
T _1 reflects interference between 2 terms of equal amplitudes 
and densities (n sl , n a2 ) = (1.8, 1.7) xlO 12 cm" 2 . The fit 
yields mobility fx = 2,800 cm 2 /Vs, and k F l = 41. Panel (b) 
shows the fit of the envelope versus T with B fixed at 12 T. 
The fit yields m* = 0.089 m (free mass). With k F = 0.047 
A -1 , the inferred velocity vf — 6 xlO 5 m/s. 

likely arise from the large surfaces of the cleaved crys- 
tal that are normal to B. From the inferred n s i and p, 
we find that the conductance of each surface is G s = 
\{e 2 /h)k F l ~ 0.72 mS (or R D ~ 1.39 Ml). 

The prominence of the oscillations allows us to address 
the question whether the surface states have a Dirac dis- 
persion. In principle, one may plot the "index" field 1/B V 
versus the integers v = 1, 2 • • • , and track the intercept in 
the limit B — ► do. However, away from the quantum Hall 
effect (QHE) regime at low B, it is sometimes uncertain 
whether one should take the maxima or minima of p xx 
or p yx for the index field. To us, the most natural choice 
is the field defined by the filling factor v = N e /N$, viz. 

B v = n s <p /v, (3) 

where N e and are the total number of electrons on a 
surface and the number of flux quanta 4>q = h/e pierc- 



ing the surface (hereafter, we focus on one surface, i.e. 
one spin degree of freedom) . Equation [3] is equivalent to 
having exactly v flux quanta enclosed in the cyclotron 
orbit for an electron in the LL v, viz. Birkp£g = v<f>$, 
where £b = y/h/eB is the magnetic length. Hence, the 
chemical potential £ lies between 2 LLs when B = B v . 

To show how the line of 1/B„ vs. v relates to the 
Dirac dispersion, we employ the construction in Fig. |4] 
In Panel (a), we have sketched the steps of the quantized 
Hall conductivity a xy — (e 2 /h){v + |) versus energy E, 
where the important shift of \ arises from the v = 
LL at the Dirac point, as in graphene. The half-ovals 
drawn on the iJ-axis represent the broadened LLs. As 
the peaks occur at the step-edges of a xyi they also locate 
the maxima in dp yx /dB. Consequently, the minima in 
\dp yx /dB\ locate B„, the field at which £ lies between 
LLs (arrows). The ^-shift means that the line of 1/B V 
vs. v in Panel (b) does not pass through the origin (unlike 
the case of a quadratic dispersion). 

In Fig. HJd, we have plotted the measured 1/B V versus 
v. With just one vertical scale adjustment, we can align 
the 5 data points to the steps in Panel (a). At our highest 
B, the filling factor is v — 5. More significantly, the line 
through the data in (b) intercepts the v axis at v = -0.55 
instead of 0. Hence our results are consistent with the | 
shift expected from a Dirac spectrum. 

Recently, Fisher et al. observed fractional-filling states 
in (Bi,Sb)Se3 in very intense fields (>50 T). In our ex- 
periment, evidence for fractional- filling begin to emerge 
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at much lower fields (11 T). As shown in Fig. [5ji, an 
intriguing array of sharp maxima in dp yx /dB (arrows) 
is apparent for 4< v <6. The peaks become weaker as 
T increases from 0.35 K, but some are still resolved at 
10 K. Interestingly, the peak positions align well with 
fractional values of v, as shown in Fig. [S]d. 



As in the fractional QHE regime in GaAs and 
graphene, we interpret the peaks as evidence for many- 
body states that are stabilized at fractional v. However, 
there are several puzzling features (when compared with 
the standard FQHE phenomenology). First, in Fig. [5ji, 
we find that the maxima of \dp yx /dB\ locate the frac- 
tional values of v whereas B v for integer v is fixed by 
the minima. We do not have an explanation for this in- 
version. Second, the fractional-filling peaks are observed 
at fairly large v, whereas in GaAs, they are difficult to 
resolve for v >2, despite the much higher fi s in GaAs. 
Finally, unlike in GaAs, the fractions corresponding to ^ 
are more prominent here (particularly 11/2) than those 
corresponding to | or |. In this system, we may in- 
corporate featurs such as the linear Dirac dispersion, 
strong spin-orbit interaction, and strong suppression of 
2kp scattering to understand the states at fractional fill- 
ing. 



FIG. 5: (Color online) Panel (a): Expanded view of dp yx /dB 
for 4< v <6 showing sharp maxima at non-integer values of v, 
at T from 0.35-10 K. The arrows locate the more prominent 
peaks. In Panel (b), the peak positions plotted against v align 
well with fractional values of v (arrows mark the values v — 
13/3, 9/2, - 17/3). 
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